Effects of Dynamical Quarks on the Stability of Heavy Quarkonia in Quark-Gluon 
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To study the stability of heavy quarkonia in quark-gluon plasma, we use a color-singlet Q-Q 
potential determined previously to be a well-defined linear combination of the free energy Fi 
and the internal energy Ui . Using the lattice gauge results of Kaczmarek et al., the dissociation 
temperatures of J/ip and xt> in quenched QCD are found to be 1.62T C and 1.18T C respectively, 
in good agreement with spectral function analyses. The dissociation temperature of J/ip in 
full QCD with 2 flavors is 1.42T C . Thus, the presence of dynamical quarks in full QCD lowers 
the dissociation temperature of J/ip, but J/ip remains bound up to 1.42T C . 
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1 Introduction 

Previously, Matsui and Satz suggested that the suppression of J/ip production can be used as a 
diagnostic tool for quark-gluon plasma^. The stability of heavy quarkonia has been the subject 
of many recent investigations 2_13 . Recent spectral analyses of quarkonium correlators in lattice 
gauge calculations indicate that J/ip is bound up to 1.6T C in quenched QCD 4 . Kaczmarek et al. 
have carried out lattice gauge calculations and obtained the free energy Fx and the total internal 
energy U\ for a col or-sin glet Q-Q pair in quark-gluon plasma, in both quenched QCD and full 
QCD with 2 flavors'"^. It is of interest to test whether the potential model using F\ and U% 
lead to results consistent with spectral function analyses in the same quenched approximation. 
If so, the potential model can be reliably used in full QCD to assess the effects of dynamical 
quarks on the stability of heavy quarkonia in quark-gluon plasma. 



2 The Q-Q Potential and Heavy Quarkonium Stability 

The most important physical quantity in the potential model is the Q-Q potential that acts 
between a quark Q and an antiquark Q. The potential should be obtained by proper theoretical 
considerations and must yield results consistent with spectral function analyses. 

For a Q-Q pair in a thermalized system at a fixed temperature and volume, the equation of 
motion for the QQ pair can be obtained by minimizing the grand potential with respect to the 
QQ wave function. Using such a variational principle in a schematic model, we find in Ref. 



that the proper color-singlet Q-Q potential is given by the internal energy XJ\ after subtracting 
out the quark-gluon plasma internal energy, as the quark-gluon plasma internal energy gives rise 
to the potential for the plasma constituents and not directly to the potential between Q and Q. 
Such a subtraction can be carried out by noting that the quark-gluon plasma energy density e 
is related to its entropy density a and pressure p by the first law of thermodynamics, and the 
quark-gluon plasma energy density is also related to the pressure by the equation of state. Thus 
after taking into account the response of the quark-gluon plasma to the polarizing presence of 
the Q and the Q in a hydrodynamical framework, the proper Q-Q potential, Uqq, is a linear 



combination of the free energy F\ and U\ given by 
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where a(T) = 3p(T)/e(T) can be obtained from the quark-gluon plasma equation of state. The 
potential f/l 1 ! is approximately F\ near T c and approximately 3iq/4 + £/i/4 for T > 1.5T C . 
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Figure 1: Parameters for Fi(r,T) for full QCD. 
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Figure 2: Fits to Fi(r, T) of Ref. 



for full QCD. 



Kaczmarek et al. recently obtained F\ and U\ in both quenched QCD - and full QCD 
with 2 flavors ^. In quenched QCD, F\ and U\ can be parametrized in terms of a screened 
Coulomb potential with parameters shown in Figs. 1 and 2 of Ref. [Hj. In full QCD with 2 
flavors, F\ and U\ can be simply represented by a color-Coulomb interaction at short distances 
and a completely screened, constant, p otentia l at large distances, although other alternative 
representations have also been presented ^^ l The transition between the two different spatial 
regions can be described by a radius parameter ro(T) and a diffuseness parameter d(T), 
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In searching for the coupling constant a s , we found that the value of a s centers around 0.3. It 
is convenient to fix a s to be 0.3 so that there are only three parameters for each temperature. 
The sets of parameters C, rg, and d for F\(y,T) and U\(r,T) in full QCD with 2 flavors are 
shown in Figs. 1 and 3 respectively, and the corresponding fits to the F\ and XJ\ lattice results of 



Kaczmarek et al. for full QCD with 2 flavors are shown in Figs. 2 a nd 4 . To calculate the ratio 
a{T) in Eq. (1), we use the quenche d e quation of state of Boyd et a for quenched QCD, and 
the equation of state of Karsch et aA^"for full QCD with 2 flavors. Using quark mases m c = 1.4 
GeV and = 4.3 GeV, we can calculate the binding energies of heavy quarkonia and their 
dissociation temperatures. We list in Table I the heavy quarkonium dissociation temperatures 
calculated with the Uqq potential, the F\ potential, and the U\ potential, in both quenched 
QCD and full QCD. 
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Figure 3. Parameters for Ux(r, T) for full QCD. Figure 4. Fits to U x (r, T) of Ref. for full QCD. 



With the Uqq potential, the dissociation temperatures of J/tp and \b are found to be 1.62T C 



and 1.18T C respectively in quenched QCD. On the other hand, spectral analyses in quenc hed 
QCD give the dissociation temperature of 1.62-1.707^. for J/ip® and 1.15-1.54T C for Xb 
Thus, the Uqq potential of Eq. Q gives dissociation temperatures that agree with those from 
spectral function analyses. In contrast, the U\ potential in quenched QCD gives dissociation 
temperatures that are much higher while the F\ potential gives dissociation temperatures that 
are lower than the corresponding temperatures from spectral function analyses. Therefore Uqq 
is the proper Q-Q potential for studying the stability of heavy quarkonia in quark-gluon plasma. 



Table I. Dissociation temperatures obtained from different analyses. 





Quenched QCD 


Full QCD (2 flavors) 


States 


Spectral Anal. 


QQ 


F x 


Ux 


QQ 


F x 


Ux 




1.62-1.7023 


1.62 T c 


1.40 T c 


2.60 T c 


1.42 T c 


1.21 T c 


2.22 T c 


Xc 


below l.lTj 


unbound 


unbound 


1.18 T c 


1.05 T c 


unbound 


1.17T C 






unbound 


unbound 


1.23 T c 


unbound 


unbound 


1.11 T c 






4.1 T c 


3.5 T c 


~ 5.0 T c 


3.30 T c 


2.90 T c 


4.18 T c 


Xb 


1.15-1.542? 


1.18 T c 


1.10 T c 


1.73 T c 


1.22 T c 


1.07 T c 


1.61 T c 


TW 6 




1.38 T c 


1.19 T c 


2.28 T c 


1.18 T c 


1.06 T c 


1.47 T c 



tRcf?" 



^ReP, "Ref 111 " 



In full QCD with 2 flavors, the dissociation temperature of J/ip calculated with the J7 
potential is 1.42T C , as listed in Table I. Thus, the presence of dynamical quarks in full Q 



lowers the dissociation temperature of J/ip because of the additional screening, but J/ip remains 
bound up to 1.42T C . 



3 Discussions and Conclusions 

The stability of quarkonium depends on the potential that acts between the qu ark Q and an 
antiquark O. P revious work in the potential model uses the free energy i^EEEa or the internal 
energy Ui " obtained from lattice gauge calculations as the Q-Q potential. 

We find in Ref. [ that the proper Q-Q potential, Uqq, is actually a well-defined linear 
combination of Fi and U\ with coefficients that depend on the equation of state. We test such 
a Q-Q potential by evaluating the dissociation temperatures for heavy quarkonia. Within the 
same quenched approximation, the UqL potential give results that are consistent with those 

from spectral function analyses. We can therefore use the Uqq potential in full QCD to study 
the effects of dynamical quarks on the stability of heavy quarkonia. 

We find that for J/ip, the dissociation temperature is changed from 1.62T C for quenched 
QCD to 1.42T C for full QCD with 2 flavors. Thus, the presence of dynamical quarks in full QCD 
lowers the dissociation temperature of J/ip because of the additional screening due to dynamical 
quarks. The change in stability is nonetheless quite small, as J/ip remains bound up to 1.42T C 
even in the presence of dynamical quarks. 

In conclusion, we have shown that the potential model is consistent with the spectral function 
analysis in quenched QCD, and we have found that the effects of dynamical quarks lowers the 
J/ip dissociation temperature to 1.42T C . 
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